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ABSTRACT 

Context. Protoplanetary disks generally exhibit strong line emission from the CO fundamental v=l~0 ro- vibrational band around 
4J fim. The lines are usually interpreted as being formed in the Keplerian disk, as opposed to other kinematic components of the 
young stellar system. 

Aims. This paper investigates a set of disks that show CO emission line profiles characterized by a single, narrow peak and a broad 
base extending to > 50 km s"' , not readily explained by just Keplerian motions of gas in the inner disk. 

Methods. High resolution {R= 10^) M-band spectroscopy has been obtained using CRIRES at the Very Large Telescope in order to 
fully resolve fundamental ro-vibrational CO emission line profiles around 4.7 yum. 

Results. Line profiles with a narrow peak and broad wings are found for 8 disks among a sample of ~50 disks around T Tauri stars 
with CO emission. The lines are very symmetric, have high line/continuum ratios and have central velocity shifts of < 5 kms"' 
relative to the stellar radial velocity. The disks in this subsample are accreting onto their central stars at high rates relative to the parent 
sample. All 8 disks show CO emission lines from the v = 2 vibrational state and 4/8 disks show emission up to v = 4. Excitation 
analyses of the integrated line fluxes reveal a significant difference between typical rotational (~300-800 K) and vibrational (~1700 
K) temperatures, suggesting that the lines are excited, at least in part, by UV-fluorescence. For at least one source, the narrow and 
broad components show different excitation temperatures, but generally the two component fits have similar central velocities and 
temperature. Analysis of their spatial distribution shows that the lines are formed within a few AU of the central star. 
Conclusions. It is concluded that these broad centrally peaked line profiles are inconsistent with the double peaked profiles expected 
from just an inclined disk in Keplerian rotation. Models in which the low velocity emission arises from large disk radii are excluded 
based on the small spatial distribution. Alternative non-Keplerian line formation mechanisms are discussed, including thermally and 
magnetically launched winds and funnel flows. The most likely interpretation is that the broad-based centrally peaked line profiles 
originate from a combination of emission from the inner part (< a few AU) of a circumstellar disk, perhaps with enhanced turbulence, 
and a slow moving disk wind, launched by either EUV emission or soft X-rays. 

Key words. Protoplanetary disks - Line: profiles - Stars: low-mass - Planets and satellites: formation - Accretion, accretion disks 

1 . Introduction the lifetime of the gas in the inner part of the disk and thus its 

ability to form giant gaseous planets. 

It is generally thought that plan ets form i n the inner regions 

of protoplanetary disks (<10 AU iLissaueii IT993I) . Information The bulk of the gas mass m protoplanetary disks is in H2 

on the physical structure, gas dynamics and chemical com- ^ut this molecule is difficult to observe since its rotational 

position of the planet-forming region is essential to constrain quadrupole transitions from low-energy levels are intrinsically 

models of planet formation. Processes like planet migration, ^^^k and he in wavelength ranges witlino or poor atmospheric 

^ which can change the orbits of newly formed planets, de- transmission (e.g., Carmona^ Ml). In contrast, the next 

■ pend sensit ively on the pre sence of gas in the disk (e.g., most abundant molecule, CO, has ro-vibrational lines which can 

■ ■ ' [&,T997l'Klev etaUlIOp. The planetary mass distributions readily detected from the ground. This makes CO an opti- 

(Ida & Lin, 2004) and planetary orbits ( Kominami & Ida, 2002; "^^1 t^^'^^'" *e characteristics oTthe warm gas in the mner 

rTrillingetal.L l2062l) resulting from planet formation models can I'egio"^ of disks (see|NajitaetaLl|2007l for overview). CO over- 

eventually be t ested against rec e nt obse rvations of exo-planetary ^0"^ emission ( Av = 2) was detected forjhe first time in low and 

systems (e.g., iMordasini et al.L l2009illbh . Observations of line ^igh mass young stellar objects by Thompson 09851) and was 

emission from disks at high spectral and spatial resolution are attributed to circumstellar disks by £ari (1989). 

needed to provide the initial conditions for these models. Gas- Overtone emission lines at 2.3 /zm from disks around T 

phase tracers of the disk surface can also b e used to probe photo- Tauri stars, when present, have been fitted with double peaked 

evapo ration processes jGorti et all 120091 ; iGorti & Hollenbacfl line profiles with a FWHM of around 100 km s which sug- 

120091) ■ More generally, these observations provide constraints on gests an origin in the innermost part (0.05-0. 3 AU) of the disk 

under the assumption of Keplerian rotation (ICarr et al 

* This work is based on observations collected at the European Chandler et al., 1993; Naiita et al., 1996). The relative intensities 

Southern Observatory Very Large Telescope under program ID 179.C- of the ro-vibrational lines can be used to determine characteristic 

0151. CO excitation temperatures (rotational and vibrational), which. 
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in turn, provide constraints on the kinetic temperatures and den- 
sities in the line-forming region. The overtone data resuhed in 
temperature estimates in the 1500-4000 K range and with den- 
sities of >10i° cm-3 dChandler et all 1 19931: iNaiita et all l2000l 
l2007h . 

Fundamental CO v = 1 - emission at 4.7 um has been ob- 

I ( 

served both from disks around T Tauri stars (e.g., Naiita et al., 
2003; Rettig et al., 2004; Salyk et al., 2007; Pontoppidan et al.,, 
2008; Salvketal., 2009) a nd around He r big A e/ Be stars (e.g 
Brittain et aL ,20 03; Blake & Booged l2004l iBrittain et all 
20071 120091; Ivan der Plas et al.. 200^ The fundamental ro- 
vibrational lines are excited at lower temperatures (1000-1500 
K) than the overtone lines (Naiita et al., 2007). Fundamental CO 
emission lines are usually fitted with double peaked or narrow 
single-peaked line profiles that can be described by a Keplerian 
model. Single-peaked line profiles with broad wing s (up to ICj O 
km s ') in T Tauri disks have been seen by Naiita e t al.l (l2003h . 
who also discussed their origin. Three possible formation sce- 
narios were mentioned: disk winds, funnel flows or gas in the 
rotating disk. The latter option was favored, and the lack of a 
double peak was ascribed to the relatively low (R=A/AA=25QQ0) 
spectral resolving power of the data. Alternative explanations 
such as an origin in disk winds and funnel flows were ruled out 
mainly because of the lack of asymmetry in the line profiles. 

The high resolution (R = 10^) spectrometer CRIRES 
(CRyogenic InfraRed Echelle Spectrograph) fed by the MACAO 
(Multi - Application Curvature Adaptive Optics) adaptive optics 
system on the Very Large Telescope offers the opportunity to 
observe molecular gas emission from T- Tauri disks with unsur- 
passed spectral and spatial resolution. A sample of ~70 disks 
was observed in the fundamental CO band around 4.7 fim as part 
of an extensive survey of molecular emission from young stellar 
objects. In total, 12 of the 70 T Tauri stars show CO emission 
lines with a broad base and a narrow central peak, from now on 
called broad-based single peaked line profiles. Eight of the 12 T 
Tauri stars are selected for detailed analysis in this paper, based 
on criteria discussed in ^3.51 Because the lines remain single 
peaked even when observed at 4 times higher spectral resolu- 
tion than previous observations, the lack of a double peak can 
no lon ger be explained by the limited resolution in Naiita et al. 
( 2003 ). Hence the modeled double peaked lines in Najita et al., 
(l2003h are not a plausible explanation for the centrally peaked 
line profiles in T Tauri disks. The aim of this paper is to classify 
these broad centrally peaked lines and to constrain their origin. 
Since many of these sources have high line to continuum ratios 
and are p rime targets to search for molecules other than CO (e.g. 
ISalvket al..2008). a better understanding of these sources is also 
warranted from the perspective of disk chemistry studies. 

The observations and sample are presented in ^ In section 
|3] the line profiles are modeled using a Keplerian disk model 
where it is concluded that a model with a standard power-law 
temperature structure does not provide a good fit to the broad- 
based single peaked line profiles. The profiles are subsequently 
inverted to determine what temperature distribution would be 
consistent with the spectra. The origin of the emission is then 
further constrained in ^by extracting radial velocity shifts be- 
tween the gas and the star, determining rotational and vibrational 
temperatures and investigating the extent of the emission. In ^ 
the results are discussed and they are summarized in ^ 

2. Observations and sample 

A sample of 70 disks around low-mass pre-main sequence 
stars was observed at high spectral resolving power (A/AA - 



10^ or 3 km s ') with CRIRES mounted on UTl at the Very 
Large Telescope (VLT) of the European Southern Observatory, 
Paranal, Chile. The CRIRES instrument (Kiiufletal., 2004) is 
fed by an adaptive optics system (MACAO, Paufigu e et all 
12004 ). resulting in a typical spatial resolution of ~ 160-200 milli- 
arcsec along the slit. CRIRES has 4 detectors that each cover 
about 0.02 - 0.03 jum with gaps of about 0.006 yum at 4.7 jum. 
Staggered pairs of settings shifted in wavelength are observed to 
cover the detector gaps and produce continuous spectra. 

The observations were taken during a period from April 21 
2007 to January 3 2009. The parent sample is a broad selec- 
tion of low-mass young stellar objects, consisting mostly of T 
Tauri stars and a few Herbig Ae stars. The full data set will 
be published in a future study (Brown et al. in prep.). Of the 
70 sources, about 50 disks around T Tauri stars show clear CO 
emission lines. This paper focuses on a subsample of 8 of the 12 
objects that show broad-based single peaked CO ro-vibrational 
line profiles. Their names and characteristic parameters are pre- 
sented in Table[T] and their selection is justified in §3.2 and §4. 
Inclinations are unknown for the majority of the sources. Several 
of the sources in the sample are binaries, with the primary and 
secondary defined as A and B, respectively. The specific def- 
inition for each source's primary and secondary is taken from 
the literature, see Table [T] For S CrA, both A and B compo- 
nents show broad-based centrally peaked line profiles. For AS 
205 and V V CrA, only the A component has such single-peaked 
emission profiles. The B components show CO absorption (see 
ISmifli et alj|2009l for the case of VV CrA B). Thus, there is no 
obvious trend of the presence of these profiles with binarity of 
the system. 

Ground-based M-band spectroscopy is usually dominated by 
strong sky emission lines (including CO itself) superimposed on 
a thermal continuum. Nodding with a throw of 10" was per- 
formed to correct for both the sky emission and the thermal 
continuum. In addition jittering with a random offset within a 
radius of Of.'S was used to decrease systematics and to remove 
bad pixels. The weather conditions during the observations were 
usually good with optical seeing varying typically between 0'.'5 
- 1'.'5. Even at times when the seeing was higher most of the 
light still passed through the 0'.'2 slit because of the adaptive op- 
tics system. Observations of standard stars were done close in 
time to each science target with airmass differences of typically 
0.05-0. 1 in order to correct for telluric features. 

The observations typically covered CO V - 1-0: R(8)-P(32) 
lines, where the notation indicates that most lines between R(8) 
and P(32) are observed, including, for example, the R(0) and 
P(l) lines. Lines from vibrationally excited levels as well as iso- 
topologues are also included, specifically: CO v = 2 - 1: R(17)- 
P(26); CO V = 3 - 2: R(26)-P (21); CO v = 4 - 3: R(37)-P(15); 
'^CO: R(24)-P(23); C'^O: R (17)-P(27); and C'^O: R(26)- 
P(22). 

The wavelengths and dates of the observations can be found 
in Table |2l Each wavelength setting is centered on the third, in 
order of increasing wavelength, of the four detectors. As an ex- 
ample, the setting at 4.730 /vm corresponds to a spectral range of 
4.660 - 4.769 yum. Observations were sometimes taken during 
different seasons to shift the telluric lines relative to the source 
spectrum due to the reflex motion of the Earth. Combination of 
the two data sets then allows the reconstruction of the complete 
line profile, whereas a single epoch would have a gap due to the 
presence of a saturated telluric CO line. For 6 of the sources pre- 
sented here, no evidence has been found that the line fluxes and 
profiles changed within this 1-1.5 yr period. For VW Cha and 
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Table 1. List of sources with single-peaked line profiles analyzed in this paper. 



Source" 


Q'(J2000) 


^(J2000) 


Spectral Type 


Distance (pc) 


Flux [Jy] 


Teff [K] 




Ref." 


AS 205 A (N:i;_ 


16 11 31.4 


-18 38 24.5 


K5 


125 


4.0 


4250 


4.0 


1,4, 8, 10 


DRTau 


04 47 06.2 


+ 16 58 42.9 


K7 


140 


1.3 


4060 


1.1 


4, 5, 11 


RULup 


15 56 42.3 


-37 49 15.5 


K7 


140 


1.1 


4000 


1.3 


1,6, 12, 15 


S CrA A (N) 


19 01 08.6 


-36 57 20.0 


K3 


130 


2.2 


4800 


2.3, 


2, 3, 17 


S CrA B (S) 


19 01 08.6 


-36 57 20.0 


MO 


130 


0.8 


3800 


0.8 


2, 3, 17 


VV CrA A (S) 


19 03 06.7 


-37 12 49.7 


K7 


130 




4000 


0.3 


3, 13, 16 


VW Cha 


11 08 01.8 


-77 42 28.8 


K5 


178 


0.7 


4350 


2.9 


7, 9, 14 


VZ Cha 


11 09 23.8 


-76 23 20.8 


K6 


178 


0.4 


4200 


0.5 


7, 9, 14 



" Three of these sources are binaries and the separations between 
their A and B components are: A S 205: 1'.'3, S CrA: 1'.'3 and VV CrA: 
V:9 <Reipurth & Zinneckerl.ll993h. 

* Re ference s. - ( 1) 'Kessler-Silacci et al.' ('2006'); (2) 'Pr ato et alj 
(l2003h: G) iTakam ietal. (2003); (4) Salvketal. (200^); (5) 
Muzerolle et al. (2003); (6) Gtinther & Schmitt ( 200 ^, (7) 
Luhmanetal. (2008), (8) Evanset al. (2009), (9) Whitt et etalj 
ll997l); (10) Andrews et aU feOOgI) ; (11) Ricci et al. (20l3); 
(12) iGrag-Ve lazguez & Ray, ( l2005h; (13) ,Koresko et al.. (19971) ; 
(14) Natta et al. (2000); (15) IStempels & Piskunovl ( |2003[) (16) 
lAppenzeller et al.i (il986) and (17) Peterson et al. (subm.). 

N or S indicate if the source is the northern (N) or southern (S) of a 
binary system. 



Table 2. Journal of observations. 



Source 


Obs. time 


Settings (yum^ 


Standard star 


Spectral type 










standard star 


AS 205 A 


Apr 07 


4.760, 4.662, 4.676, 4.773 


BS 4757 


AO 




Aug 07 


4.730 


BS 5812 


B2.5 




Apr 08 


4.730 


BS 6084 


Bl 




Aug 09 


5100, 5115 


BS 5984 


B0.5 


DRTau 


Oct 07 


4.716, 4.730, 4.833,4.868 


BS3117, BS 838 


B3, B8 




Dec 08 


4.716, 4.946 


BS 1791 


B7 


RULup 


Apr 07 


4.716, 4.730, 4.833, 4.929 


BS 5883 


B9 




Apr 08 


4.730 


BS 6084 


Bl 


SCrA 


Apr 07 


4.730, 4.716 


BS 6084 


Bl 




Aug 07 


4.730 


BS 7235, BS 7920 


AO, A9 




Aug 08 


4.868, 4.946 


BS 6084, BS 7236 


Bl, B9 


VV CrA A 


Apr 07 


4.716, 4.730, 4.840 


BS 7362 


A4 




Aug 07 


4.770, 4.779 


BS 7236 


B9 




Aug 08 


4.946 


BS 7236 


B9 


VWCha 


Dec 08 


4.716, 4.800, 4.820, 4.946 


BS 5571, HR4467 


B2, B9 


VZ Cha 


Dec 08 


4.716, 4.800, 4.820, 4.946 


BS 5571, HR4467 


B2, B9 



" The reference wavelength of a given setting is centered on the third 
detector. 



VZ Cha this information is not available since they were only 
observed once. 



2.1. Data reduction 

The spectra were reduced using standard m ethods for infrared 
spectroscopy (see Pontoppidan et al.L 12008 1 for details). Dome 
flats were used for each setting to correct the images for pixel-to- 
pixel sensitivity differences in the detectors. A linearity correc- 
tion was applied to the frames using the parametrization given 
in the CRIRES documentation. The nod pairs were then dif- 
ferenced to subtract the background. The data stack was co- 
added after correction for field distortion. Extraction of the 2- 
dimensional images to 1 dimensional spectra was done after the 
subtraction and combination of the 2-D frames tracing the source 



in the 2-D spectra. Additional spectra of well separated binaries 
were extracted as well during this procedure. 

The spectra were subsequently wavelength calibrated by fit- 
ting the standard star atmospheric lines to an atmospheric model 
generated using the Reference Forward Model (RFM) code. This 
model is a line-by-line rad iative transfer model developed at 
Oxford University based on 'Clough et al. ' ('198213 The typical 
velocity accuracy is about 0.1-1 km s"', but can vary between 
settings, depending on the density of telluric lines. The last step 
in the reduction process was to correct for the strong telluric 
absorption lines in the spectra, by dividing the science spectra 
with the spectra taken of early-type photospheric standard stars. 
These hot stars have a strong hydrogen Pfund fi (7-5) absorption 
line at 4.654 jim. Division by the standard star therefore intro- 



' |http : //www ■ atm . ox . ac . uk/RFM/] 
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4.630 4.640 4.650 4.660 



Wavelength [/^m] 

Fig. 1. Part of the VLT-CRIRES spectrum of AS 205 A. The strong '^CO Hnes are labelled with their /?-band lower 7-transition 
number and the other marked lines show detections of '^CO v = 2-1, v = 3-2, '^^CO, and the Pfund-/? line. The gaps in the spectrum 
are excised telluric features. For other parts of the AS 205 A spectrum, these gaps have been filled in by observations at different 
times of the year (see Fig.|2ll. 



duces an artificial contribution to the emission line in the sci- 
ence spectra. However, Pfund fi is detected in emission towards 
most of the sources prior to telluric correction. Early-type photo- 
spheric standard stars are otherwise rather featureless throughout 
the spectral band. Flux calibration was carried out by scaling to 
photometry from the Spitzer IRAC band 2 (see Table[T]). 

In Fig.[T]part of the reduced spectrum for the northern of the 
binary components of AS 205 is presented. This spectrum shows 
highly resolved '^CO line profiles for the v = 1 - transitions 
R(8)-R(0) and clear- detections of '^CO v= 1 - 0, CO v = 2 - 1 
and V = 3 - 2 lines and the Pfund fi line. An overview of parts of 
the spectra of the 8 single peaked sources is presented in Fig.|2] 



3. Line profiles 

3.1. ^^CO line profiles 

The spectrally resolved fundamental CO emission lines from the 
sample of T Tauri stars observed with CRIRES show a broad 
variety of profiles. Often, the line profiles can be explained by 
Keplerian rotation, absorption by the disk or absorption by a 
foreground cloud (Brown et al., in prep). Fig.[3]shows a selection 
of emission line profiles of three basic shapes: A) narrow single- 
peaked, matched by a single Gaussian, B) double-peaked line 
profiles and C) single-peaked with broad wings. Lines that show 
absorption are not included due to the difficulty in determining 
the profile close to line center. 



The CO line profiles of TW Hya fall in category A and 
likel y arise from g a s in a face-on disk in Keplerian rotation (e.g. 
Pon toppidan et al. , |2008|) . Prototypical for category C are the 
line profiles of the source AS 205 A, which are by inspection 
qualitatively different from the double peaked profiles of sources 
like AA Tau and VV Ser (category B) because of their single 
narrow peak relative to a broad base. While categories A and B 
appear to be well understood (see below), the question is what 
the origin of the broad centrally peaked line profiles of category 
Cis. 

3.2. Kepierian disk modei 

To investigate to what extent the broad centrally peaked (cat- 
egory C) line profiles can be described by emission origi- 
nating from gas circulating in the inner parts of the disk, a 
simple d isk mode l of gas in Keplerian rotation adopted from 
IPontoppidan et al.l (i2008.) is used. 

3.2.1 . Description of tine model 

The model includes a flat disk with constant radial surface den- 
sity. The motions of the gas in the disk are described by Kepler's 
law. The disk is divided into emitting rings with radial sizes that 
are increasing logarithmically with radius. Each ring has a spe- 
cific temperature. Local thermal equilibrium (LTE) is used as 
an approximation for describing the excitation of the emitting 
molecules, which means that the relative populations in a ro- 
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DR Tau 






RU Lup 





S CrA A 




S CrA B 








2.5 
2.0 
1.5 
1.0 






Wavelength [/xm] 

Fig. 2. Partial spectra for the entire sample of sources with single peaked '^CO lines. The HI Pfund /3 line at 4.655 /im is included. 
The spikes seen in spectra of S CrA A, B and VV CrA A are narrow (self) absorption lines visible in the low J lines. 



5 



J.Bast et al.: Single peaked CO emission line profiles from inner regions of disks 




-40 -20 20 40 -100 -50 50 100 -50 50 



Velocity [km s'^] 

Fig. 3. Modeled lines, plotted in red, are fitted to observed spectral lines in black of TW Hya, VV Ser and AS 205 A. Good fits could 
be found for VV Ser and TW Hya using a Keplerian model with a power-law temperature profile. No good fits could be found for 
AS 205 A. The optimal fit to the total line profile (red) using extended emission out to 10 AU, to the narrow inner part (green) and 
to the outer broader wings (blue) are presented in the right figure. See text for model parameters. 



vibrational level can be derived using the Boltzmann distribution 
at the given temperature in each emitting ring. 

In the standard formulation, a power law is used to describe 
the temperature gradient throughout the disk, T - ToiRIRoY" , 
where T is the gas temperature, R is the radius and Tq is the 
temperature at the inner radius Rq. The other parameters are: the 
mass of the star ; the inclination of the disk /; the outer radius 
of the disk Rmax, and a line broadening parameter AVmod which 
combines the instrumental broadening and any turbulent broad- 
ening. Several parameters are degenerate, including the stellar 
mass and disk inclination. Thus, this model is not used to derive 
the best estimates of the characteristic parameters of the sources. 
The aim is instead to explore the parameter space to see whether 
good fits can be found for the category C line profiles using ba- 
sic Keplerian physics. The local line profile is convolved with 
a Gaussian turbulent line FWHM broadening (AVturb) which de- 
scribes the local turbulence of the gas in the disk. The instrumen- 
tal broadening (AVinstr) can also be represented by a Gaussian, 
which has a FWHM value of 3 km s^' for CRIRES. The lowest 
value for AVmod = (A^t?irb ^^Lti-)'^^ therefore set to 3 km 
s"' and any additional broadening represents an increase in the 
value of AVtuib. For temperatures up to 1000 K, typical of the 
molecular gas in the inner disk, the thermal broadening of the 
CO lines of up to 1.3 km s ' is negligible. 

3.2.2. Fits with a Keplerian model with a power-law 
temperature profile 

The sources TW Hya, VV Ser and AS 205 A were chosen as 
illustration since their fundamental CO emission lines represent 
the three main different types of line profiles in our large 
sample. Each of these sources and their fits are presented here 
and shown in Fig. [3] 

Single peaked narrow line profile - TW Hya: A good fit 
was found to the narrow single peaked line profile of the TW 
Hya disk. This fit was achieved by using a low inclination angle 
of / = 4.3°, a = 0.4, a stellar mass of 0.7 Mq, /?o = 0.1 AU, Tq = 
1100 K, /?n,ax - 1.5 AU and AVmod =3 km s"', see the left plot 
in Fi g. [3] The ado pted st e llar rn ass and disk inclination are taken 
from iPontoppida n et al.l (l2008h . For this source, CO emission 
from the Keplerian disk model described above, convolved to 
the 3 km s ' instrumental resolution of CRIRES, provides a 
good fit to all lines. Additional line broadening from turbulence 
is not needed to explain the line profiles. This low turbulence 



is consistent with the low turbulent velocity of ~0.1 km s 
inferred for the outer disk of TW Hya by lOi et alJ (l2006l) . 
although that value refers to much cooler gas at large radii in the 
disk. 

Double peaked line profile - VV Ser: VV Ser is an 

example of a disk with double peaked line profiles extending 
out to velocities of + 40 km s"'. VV Ser is a Herbig Ae/Be 
star of spectral type A2-B6 with a mass of 2.6 + 0.2 Mq and 
an inclination of 65-75" (Ponto ppidan et all l2007h . A good fit 
between the model and the data is shown in the center panel in 
Fig[3] The estimated parameters are M* = 2.6 M©, / - 70°, Tq - 
3500 K, ff = 0.45, R() = 0.08 AU, R^^^ = 1 1 AU and AV^od =3 
km s"' . However in this case the parameter space is rather large 
since smaller inclination angles will give equally good fits if a 
higher stellar mass is used. The main point here is that a good 
fit can readily be obtained with reasonable parameters for a 
Keplerian model. 

Single peaked broad-based line profile - AS 205 A: No 

equally good fit could be found for the emission lines from AS 
205 A. An approximate fit is shown in red in the right plot of 
Fig. [3] with a model using an inclination angle of / = 22°, Rq - 
0.04 AU, /?max = 5 AU, AVmod =3 km s"', Tq = 1100 K, a = 
0.3, and a mass of 1.0 M©. The ma ss is the same and the inclina- 
tion angle is close to that given bv lAndrews et al.l (l2009l 1.0 Mq 
and 25°, respe c tively ). As for TW Hya, the millimeter data of 
[Andrews et al.l (l2009l) refer to much larger radii (typically >50 
AU) than the CRIRES observations that the model is fitted to. A 
lower inclination angle of 10° gives a better fit to the narrow cen- 
tral part of the line but results in poor fitting of the outer wings of 
the profile, as presented by the green profile in Fig. |3] The blue 
line represents the best fit to the line wings, which is achieved by 
increasing the inclination angle to 27° and taking a to be 0.36. 

One option to improve the fit to the AS 205 A data could 
be to increase the turbulent broadening so that the central dip is 
filled in. Indeed, increasing AV,nod to 8 km s"', removes the cen- 
tral double peak (see Fig.|4]l. This fit is unable to simultaneously 
match both the narrow peak and the broad base of the profile, 
however. 

In summary, the fits to TW Hya and VV Ser show that both 
narrow single peaked and double peaked line profiles can be re- 
produced with a simple Keplerian model with reasonable param- 
eters. However a single peaked line with broad wings cannot be 
well explained with this type of standard Keplerian model as 
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Fig. 4. Fit to the broad-based single peaked line of AS 205 A (black) using a standard Keplerian model with a power-law temperature 
profile (red) including an enhanced line broadening parameter (AVmod) of 8 km s '. The other parameters are / = 22°, Rq - 0.04 
AU, /^max = 5 AU, To = 1 100 K, a = 0.3 and a mass of 1.0 Mg (as in gtSJi. 



long as the temperature gradient is described by a continuous 
power-law. 

3.2.3. Keplerian disk models with a non-standard 
temperature profile 

Since disk models with a standard temperature and density 
power-law do not fit the data, one complementary approach is 
to investigate what physical distribution would be needed to re- 
produce the line profiles within a Keplerian model. For this, an 
iterative approach is used in which the observed line profile is 
'inverted' to determine the temperature distribution that would 
be consistent with the data. Specifically, the disk is divided into 
rings and the temperature of each ring is adjusted to match the 
line profile. The left panels of Fig.|5]present the model compared 
with the velocity profile of the AS 205 CO v = 1-0 P(8) Une 
whereas the right plots show the corresponding excitation tem- 
perature profiles with radius. For LTE excitation and optically 
thin emission, the latter can also be viewed as an intensity pro- 
file with radius. The upper plots use the same parameters as in 
gTZ2l for AS 205 A but with a higher stellar mass of 1.4 Mq, 
a higher inclination angle of /=25°and including extended emis- 
sion out to 30 AU. The lower plots show a model using the same 
parameters as in ^3.2.2l but with an extended emission out to 10 
AU and increasing the AVmod to 7 km s"'. Both of these models 
fit the line profiles well, illustrating that a Keplerian disk model 
with an unusual temperature distribution can explain the data. 
The emission is much more extended, however, than in the case 
of the standard disk models with a power-law temperature pro- 
file: out to 30 AU in the case of no turbulence or out to 10 AU 
with increasing AVmod- These models will be further tested in 
95.1! using the constraints on the spatial extent of the observed 
emission found in ^4.31 



3.3. Line profile parameter 

A parameter is defined to quantify the difference in line profiles 
between the broad single peaked sources (category C) and the 
other sources in the sample. This so-called line profile parameter 
Pio describes the degree to which line profiles have a broad base 
relative to their peak. The line profile parameter Pio is there- 
fore defined as the full width (AVio) of the line at 10% of its 
height divided by the full width at 90% (AVgo) of its height, Pw 
= AVio/AVgo. The broad-based single peaked lines (category C) 
have a higher value of the line profile parameter relative to the 
double and narrow single peaked lines, see Fig. |6] 



A summary of the line profile parameter values for 31 se- 
lected T Tauri stars from the total sample of ~50 T Tauri stars 
with CO emission is presented in Fig.Qand Table|3] These 31 T 
Tauri stars were selected because their CO emission line profiles 
have high S /N and are not contaminated by strong telluric or ab- 
sorption lines. Figure|7]shows that the double peaked (turquoise 
diamonds) and narrow single peaked (green stars) sources all 
have a line profile parameter of <6. For reference, a Gaussian 
line profile has a Pio-value of 4.7. The selected sample of 8 
broad centrally peaked sources (red triangles) all have a Piq- 
value of >6. An overview of the normalized and continuum- 
subtracted '^CO P(8) line profiles of this subset is presented in 
Fig. |8] These lines are typically symmetric around line center 
and have a narrow top and a very broad base that extends out to 
+ 100 km s"' in several cases. Four other sources, DoAr24E A, 
DoAr44, EX Lup and HD135344B, also have higher (>6) Pio- 
values but are not selected for our sample (see 33.5! for source 
selection criteria). 

3.4. Modei iine profiie parameters 

To derive the maximum line profile parameter consistent with a 
standard Keplerian model with a power-law temperature profile, 
a grid of Keplerian line profiles was constructed. Three param- 
eters in the model were set to fixed values that give strongly 
peaked line profiles; Tq - 1 100 K, a - 0.30 and AVmod - 3 km 
s"'. The following parameters were varied between - 0.5 - 
1.0 Mo, / = 10 - 80°, Ro = 0.04 - 1.0 AU and R^^^ = 5-100 AU 
and the f lo- value was calculated for each produced line pro- 
file. An overview of the Pio-values versus inclination is shown 
in Fig. |9] It is found that a typical standard Keplerian model 
with a power-law temperature profile can achieve a maximum 
Pio-value of about 6, with the highest values found at low in- 
clination. This conclusion is unchanged if AVmod is increased to 
values as large as ~10 km s"'. The upper limit of f lo ~ 6 that 
the Keplerian model can reproduce is presented as a red line in 
Fig- 13 Figure|9]also indicates that the inclination of the source, 
invo ked to explain the differ ent widths of Herbig Ae disk pro- 
files dBlake&Boogertl 12004) . should not affect the value of Pio 
significantly except at very low inclinations, since it enters both 
AVio and AVgo. 

Note that the Keplerian models always require a low incli- 
nation angle of 20 - 30° to obtain high values of the line pro- 
file parameter f lo, regardless of whether or not a large turbu- 
lent broadening is included. If inclination angles within 5° of the 
range ; = 20 - 30° are considered necessary for the broad single 
peaked line profiles, then 6 + 3 sources of the total sample of 
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Fig. 5. Left panel: Keplerian model (red) using an iterative inversion method to fit the velocity profile of the AS 205 A CO v = 1 - 
P(8) line (black). Right panel: The excitation temperature as a function of radius needed to create the modeled lines. The models in 
the upper panels include extended emission out to 30 AU, a stellar mass of 1 .4 M©, an inclination of 25° and AVmod-^ km s"' . The 
models in the lower panels include emission out to 10 AU, a stellar mass of 1.0 Mq, an inclination of 22° and an increased Ayn,od=7 
km s"'. 



50 T Tauri stars would have broad, single-peaked line profiles, 
which is consistent with our subsample of 8 sources. Millimeter 
interferometry data for AS 205 A and PR Tau y i eld incUnations 
of25 ° and37+3° respectively dAndrews L [20081: [Andrews et al.L 
l2009Hlsella etall 12009) . which are roughly consistent with the 
model requirements. 

3.5. Pm-value versus the line-to-continuum ratio and source 
selection. 

The broad-based single peaked sources are also noteworthy for 
having high line-to-continuum (L/C) ratios in many lines, in- 
cluding the r esolved CO ro- vibr ational lines discussed here and 
H2O lines in lSalvketal](l2008l) . Here, the L/C-ratio is defined 
as the ratio between the peak line flux relative to the continuum 
flux. The L/C ratios of the '^CO lines have been plotted against 
the f 10-values in Fig. (TO] Typical uncertainties of the line-to- 
continuum ratio are +0. 1 . Seven of the sources stand out clearly 
with high line profile parameters and high line-to-continuum ra- 
tios. Figure [To] is used to set the criteria (Pio > 6 and L/C > 
1.6) for the selection of the broad-based single peaked line pro- 
files (category C). The lower limit for the Pio-value (f 10 > 6) for 
the sample is set to match the maximum line profile parameter 
for a line that the standard Keplerian model with a power-law 
temperature profile can produce (§3.4). The line-to-continuum 
constraint of L/C > 1.6 is rather arbitrary, based on the observa- 
tion that the majority of the sources have a lower value. Two 
other sources, EX Lup and VV CrA A, have L/C ratios less 
than 1.6 but have such a high line profile parameter (>10) that 
they are considered to belong to the category C profile sample. 
However, EX Lup is a variable source that underwent an outburst 
in 2008. The variation of the CO line profiles from EX Lup dur- 
ing and after the outburst are discussed in Goto et al. (subm.). 
DoAr 44 and DoAr24E A are borderline cases and are not in- 
cluded here because their profiles are contaminated by absorp- 
tion, limiting the analysis. HD135344B is also excluded from 



the sample due to the combination of having both a low L/C 
ratio and a Pio-value close to 6. In addition HD135344B is a 
well-known face-o n disk whose profile has been well fitted with 
a Keplerian model dPontoppidan et aU '2008; Br own et aU[2009l : 
[Gradv et al.L[2b09[) . Line profiles emitted by disks with a low in- 
clination angle can have a line profile parameter close to 6, as 
shown in Fig. [9] These selection criteria yield in total 8 sources 
that are included in the broad-based single peaked sample; these 
sources are marked with red triangles in Fig. [lO] 

4. Characteristics for the sources with broad single 
peaked lines 

In §3, we showed that our selected sources with broad single 
peaked line profiles (category C) cannot be reproduced well with 
a standard Keplerian model with a power-law temperature pro- 
file. In the following sections, we analyze the observational char- 
acteristics of the emission to obtain further constraints on the 
origin. 

4. 1 . Line profiles of CO isotopologues and thev -2 - I CO 
lines. 

Table [4| summarizes the lines seen in the spectra of the broad 
single peaked sources, including if they are detected in absorp- 
tion, emission or both. In the richest spectrum, that of AS 205 
A, lines of '^CO, C'^O, CO v = 1 - 0, v = 2 - 1, v = 3 - 2 and 
V = 4 - 3 are all detected. All of the sources, besides VZ Cha, 
have detections of '^CO. Four sources have a C'^'O detection in 
emission. In addition every source in the sample has detections 
of CO V = 2 - 1, 7 of 8 sources of CO v = 3 - 2 and 4 of 8 
sources of CO v = 4 - 3. 

Fig. [TT] compares the stacked, normalized and continuum- 
subti-acted fine profiles for '^CO v = 1 - 0, '^CO v = 1 - 0, '^CO 
v = 2 - 1 and '-CO v = 3 - 2. All of the isotopologues and the 
higher ro-vibrational transitions also show a broad-based single 
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Fig. 6. The line profile parameter Pio=Ayio/Ay9o is presented here for a sample of selected sources. Note that this parameter 
is clearly larger for the broad-based single peaked line profiles. The f lo-value is an average of the lower 7-transition lines (up 
to P(14)). The 7-transition corresponding to each plotted profile is given in the figure. The narrow absorption line seen toward 
DoAr24E A is due to '^CO absorption in the foreground cloud. 



peaked line profile. However, the profiles of the '^CO and ^^CO 
V = 1 - lines match exactly for only one source, AS 205 A. 
For the other sources, DR Tau, RU Lup, and S CrA A, the '^CO 
lines are narrower than the '^CO v = 1 - lines. 



The width of the CO v = 1-0 lines at 10% of their height 
(see Table |6ll are in general narrower than the CO v = 2 - 1 and 

V = 3 - 2 lines by about 60-80% which can be seen in Fig.fTTl 
iNaiita et"an (l2003h also found a similar trend of broader higher 
vibrational lines in their data. The difference in width between 
the lower and higher vibrational lines in our sample may reflect a 
physical difference in the location of the emitting gas. The high 
fraction of CO v = 2 - 1 and 3-2 detections in our sample 
of broad-ba s ed sin gle peaked line sources is consistent with the 
iNajita et"an (l2003h sample, where 9 of their 12 CO sources have 
detections of v = 2 - 1 emission and 3 out of 12 sources have 

V = 3 - 2 detections. However it is difficult to categorize their 
line profiles due to the lower spectral resolution. 



4.1 .1 . Two component fits 

The line profiles shown in Fig. [8] are generally not well fit by 
a single Gaussian profile. A good example of this is the CO 
emission line profiles of S CrA B, which clearly consist of two 
components. An interesting aspect is that the narrow component 
stands out very prominently in the low /-transitions but progres- 
sively decreases in intensity with higher /-transitions relative to 
the broad component (Fig.fT2li. The simplest explanation is that 
the components arise from different locations where the narrow 
component has a lower rotational temperature than the broad 
component. This is not as clearly seen for the other 7 sources 
in the sample. However, Najita et al. (2003) see a similar phe- 
nomenon in the broad-based single peaked source GW Ori. 

The line profiles from S CrA A and GW Ori support a hy- 
pothesis in which two physical components combine to form 
broad centrally peaked line profiles. We therefore fit the line pro- 
files of all sources with two Gaussians to investigate the possi- 
bility that the broad-based centrally peaked line profiles consist 
of two physical components. In the fits, all 6 Gaussian parame- 
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Fig. 7. The line profile parameter Pio for different sources. The index numbers represent the different sources, see Table 1 . Each type 
of line profile is given by a specific symbol. Data points without a symbol are for sources without clearly defined line profiles. The 
red line gives the maximum line profile parameter achieved using a standard Keplerian disk model with a power-law temperature 
profile. 
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Fig. 8. The P(8) '^CO line is plotted for the entire sample broad-based single peaked sources, apart from VW Cha for which the 
P(10) line is presented because the P(8) line is affected by a strong telluric feature at the time of observation. The lines have been 
continuum subtracted and normalized. The asymmetry of the VV CrA A line is caused by line blending by the CO v = 3 - 2 R(5) 
line. This line overlap adds intensity on the blue side of the '^CO P(8) line. 



ters (the two widths, amplitudes and central wavelengths) are left 
free. For example, Fig.[T2]shows that the CO v = 1 - P(7) line 
from DR Tau is well-fit with two Gaussian profiles, one narrow 
and one broad. 

Two Gaussian profiles are needed for 7 of the 8 broad-based 
single peaked sources to find an optimal fit to their line profiles. 
The resulting FWHMb, FWHMa,, Vb(CO) and V/vCCO) are pre- 



sented in Table|5] These results show that the broad-based single 
peaked line profiles (category C) can be fitted using a narrow 
component with a FWHM varying between ~ 10-26 km s ' 
and a broad component with FWHM of ~ 40 - 100 km s The 
central velocities of the two components are generally the same 
within the errors. Only the line profile of VZ Cha was adequately 
fit with a single Gaussian. 
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Table 3. The peakiness parameter Pio, the accretion luminosity and the line/continuum (L/C) ratio for a selected sample of proto- 
planetary disks. 



Index number 


Source 




L/C 


Acc. lum. 
[log(Lo)] 


Ref.^ 


1 


AATau 


1.8 (0.2) 


1.5 


-1.6 


1 


2 


AS 205 A 


11.5 (1.0) 


2.1 


0.2^ 


2 


3 


CV Cha 


4.1 (0.8) 


1.3 


- 


- 


4 


CWTau 


3.0 (0.2) 


1.4 


-1.3 


2 


5 


DFTau 


4.2 (0.5) 


1.5 


-0.7 


3 


6 


DRTau 


7.7 (0.8) 


2.9 


-0.1 


2 


7 


DoAr24E A 


7.3 (2.0) 


1.5 


- 1.6 


4 


8 


DoAr 44 


7.7 (0.3) 


1.4 


- 


- 


9 


EXLup 


10.8(1.0) 


1.3 


- 


- 


10 


FNTau 


3.3 (0.5) 


1.5 


- 


- 


11 


GQ Lup 


3.5 (0.5) 


1.4 


- 


- 


12 


Haro 1-4 


3.3 (0.5) 


1.4 


- 


- 


13 


Haro 1-16 


3.4 (0.4) 


1.3 


- 


- 


14 


HD135344B 


6.9 (0.3) 


1.2 


-0.9 


5 


15 


HD 142527 


6.0 (0.5) 


1.2 


0.0 


5 


16 


IRS 48 


1.7 (0.2) 


1.2 


- 


- 


17 


IRS 51 


4.1 (0.4) 


1.2 


- 


- 


18 


LkHa 330 


5.2 (0.5) 


1.3 


- 


- 


19 


RNO 90 


4.5 (0.6) 


1.5 


- 


- 


20 


RULup 


15.1 (3.5) 


1.8 


-0.4 


3 


21 


RYLup 


3.3 (0.5) 


1.2 


- 


- 


22 


S CrA A 


13.4 (2.5) 


2.0 


- 


- 


23 


S CrAB 


18.8(2.3) 


1.8 






24 


SR9 


2.4 (0.6) 


1.1 


-1.4 


4 


25 


SR21 


2.4 (0.3) 


1.2 


<-1.9 


4 


26 


TWHya 


5.0 (0.5) 


2.3 


-1.4 


3 


27 


VSSGl 


4.7 (1.0) 


1.5 


-0.4 


4 


28 


VV CrA A 


14.0 (3.0) 


1.4 






29 


VV Ser 


1.5 (0.2) 


1.1 


1.2 


5 


30 


VW Cha 


8.1 (0.5) 


2.0 


-0.2 


6 


31 


VZCha 


7.7 (0.6) 


2.2 


- 1.1 


6 



" Values refer to an average of the lower J- transitions. Uncertainties 
are indicated in parentheses 

* References - (1) GuUbrins et al] (Il998h: (2) IValenti et al.l lll993h : 
(3) iHerczeg & Hillenbrand qOO^): (4) iNatta et al.1 ( l2006h ; (5) 
iGarcia Lopez et al.1 l l2006h ~and (6) lHartmann et al.1 ( 119981) . 
For both A and B component 



Table 4. List of CO isotopologues and the '^CO v = 2- l,v = 3- 2 and v = 4 - 3 lines. 



Source" 


'^CO 


I'CO 


CO V = 2 - 1 


CO V = 3 - 2 


CO V = 4 - 3 


c'«o 


AS 205 A 


X 


X 


X 


X 


X 


X 


DRTau 


X 


X 


X 


X 


X 


X 


RULup 


X 


X 


X 


X 






SCrAA 


X, a 


X 


X 


X 


X 


X 


SCrAB 


X, a 


X 


X 


X 






VV CrA A 


X, a 


X, a 


X 


X 


X 


a 


VW Cha 


X 


a 


X 








VZ Cha 


X 




X 


X 







" X = detection in emission, a = detection in absorption and - = no 
detection. 



The results presented in Table |5] are based on fits to the CO 
V = 1 - lines that are chosen to be as uncontaminated from line 
overlap as possible. The uncertainties in Table|5]are the standard 
deviations of the fit parameters of individually fitted lines. The 
spread within a parameter for a given source is caused by un- 
certainties in the Gaussian fits, by overlapping weaker lines that 
may affect the width and position of the broad component, by er- 
rors of around ~1.0 km s"' in the wavelength calibration (which 



can vary between different detectors, see § 12.11 ). and degenera- 
cies of the fits between the broad and the narrow component. 
The large amount of line overlap in VV CrA A makes measure- 
ments of the broad component highly unreliable so the fits are 
not presented here. 

Tableland Fig. [13] show that the relative width of the peak 
and base as represented in Pio is approximately constant with 
increasing J for most of the selected category C sources, with 
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Fig. 9. Line profile parameter Pio versus inclination for a large model grid of standard Keplerian models with a power- law temper- 
ature profile. Note that Piq does not become larger than 6 and shows no trend with inclination for / > 20°. 

Table 5. The FWHM of the broad and narrow components and their average heliocentric velocities Vb{CO) and VOvCCO) relative to 
literature values of the heliocentric velocities for the stars y/,(star). 



Source 


FWHMb 


FWHMa- 


Vb(CO) 


Va,(CO) 




Ref._^ 




km s"' 


Icm s"' 


km s"' 


km s"' 


km s"' 




AS 205 A 


61.3 (16.3) 


14.8 (1.7) 


-4.4 (4.8) 


-6.2 (0.3) 


-9.4(1.5) 


1 


DRTau 


39.8 (6.3) 


13.1 (0.9) 


26.9 (1.3) 


24.4 (0.7) 


27.6 (2.0) 


3 


RULup 


96.8 (11.0) 


24.0 (3.9) 


0.4 (7.0) 


-3.9(1.1) 


-0.9 (1.2) 


1 


SCrAA 


50.2 (6.7) 


10.7 (1.2) 


-9.2(1.4) 


-5.0 (0.7) 


0.9 (0.9) 


2 


SCrAB 


97.6(17.7) 


12.3 (0.5) 


-4.7 (3.5) 


-3.5 (0.3) 


0.9 (0.9:i^ 


2 


VV CrA A;_ 




19.3 (2.0) 




-0.9 (0.2) 






VWCha 


70.6 (3.1) 


26.6 (2.2) 


15.4(1.4) 


15.1 (1.9) 


17. 2 (2.0) 


2 


VZ Clia'' 


45.1 (2.6) 




19.1 (1.7) 




16.3 (0.6) 


2 



° References. - (l) lMeicl l l2003h : (2) iGuenther et al.l l l2007h and (3) 
lArdila et al.l ^2003) . 

* Taken to be the same as S CrA A since in the literature only a value 
for S CrA is given 

^ Broad component could not be determined due to large amount of 
line overlap. 

Optimal fits with one Gaussian curve, results presented in the broad 
component columns. 



the exceptions of S CrA A and B. For these two sources the 
f lo-value decreases with increasing J (Table |6ll and the line 
profiles become dominated by the broad wings at the higher J- 
transitions, see Fig. [12] AS 205 A and VV CrA A are excluded 
in the calculations of the line profile parameter because of their 
large amount of overlapping lines that cause difficulties in reli- 
ably measuring the 10% and 90% line widths (for AS 205 A, just 
in the higher /-transitions, see Fig.[T3ll. The lack of variation in 



f 10 towards higher 7- transitions for all sources besides S CrA 
A and B, combined with the absence of a velocity shift, suggest 
that the broad and naiTOW components may be formed in related 
physical regions. 
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Fig. 10. The line to continuum ratio relative to the line profile parameter f iq. Each type of line profile is given by a specific symbol. 
Data points without a symbol do not have a clearly defined line profile. The broad single peaked sample includes sources with either 
a line/continuum ratio > 1.6 (blue line) and Pio > 6 (red line) or solely Pm > 10. 
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Fig. 11. i^CO (red), CO v = 2 - 1 (green) and v = 3 - 2 (blue) line are plotted on top of the CO v = 1-0 lines (black) for the 
broad-based single peaked CO emitting sample of T Tauri stars. Dashed lines indicate that the line profile is constructed by stacked 
lines that all are to some degree affected by line overlap. All profles have been normalized to a peak flux of unity. 



4.1.2. Central velocities 

The central velocities of the CO gas and stellar photosphere are 
compared to look for differences indicative of a wind or outflow, 
which could explain the narrow low velocity peak. The stellar ra- 
dial velocities in Table|5]are taken from literature measurements 
of photospheric optical absorption lines. The photospheric ab- 
sorption lines for our sample of sources are generally heavily 
veiled by continuum emission related to accretion, leading to 
uncertainties in the stellar velocity measurements of the order of 
l-2kms-i. 

The differences between the molecular gas radial velocity of 
the narrow component and the stellar radial velocity are small 



with shifts of IVco - Vstarl < 5 km s-\ see Fig. [H The shifts 
are more commonly seen towards the blue rather than the red, 
indicating that the CO gas is moving towards us relative to the 
star. However, the velocity shifts are only 1- 2cr and may not be 
significant. 

The small velocity differences rule out an origin of the lines 
in a fast moving disk wind or outflow. However it is still possible 
that the lines are emitted by gas in a slow disk wind. Naiit a et aP 
(2003) also do not see any significant velocity shifts between the 
radial velocity of the star and the gas within their few km s"' 
uncertainties. 
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Fig. 12. Left: The higher 7-transition '^CO P(32) (green) is plotted on top of the lower 7-transition ^^CO P(8) (black) for S CrA B. 
The narrow component of the '^CO P(8) line decreases with increasing excitation. Right: Two Gaussian fitted lines (in red and blue) 
to the CO V = 1 - P(7) line of DR Tau (black). The green line represents the sum of the two Gaussian fits. 
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Fig. 13. Comparison of the lower J and higher J profiles. The red line represents the stacked higher 7-transitions (J - P(22)-P(32)) 
and the black line the lower 7-transitions (up to 7 = P(14)) of v = 1-0 CO. The baseline for the higher stacked 7-transitions for 
AS 205 A and V V CrA A (dashed red) is not straight because of line overlap and blend. All profiles have been normalized to a peak 
flux of unity. 



Table 6. The averaged line profile parameter Pio for high and low 7-transitions for the sample of sources with single peaked CO 
V = 1 - emission lines. AVio is the averaged width of the line at 10% of its height. 



Source 


(low /) 


(high J) 


AVio [kms-'] 


AS 205 A 


11.5 (0.5J;;, 




62.8 (3.9) 


DRTau 


7.7 (0.8) 


6.6 (0.6) 


40.2 (1.0) 


RULup 


15.1 (3.5) 


10.5 (2.3) 


109(11.1) 


SCrAA 


13.4 (2.5) 


7.9 (0.2) 


56.9 (3.4) 


SCrAB 


18.8 (2.3) 


12.5 (1.0) 


102.6 (6.9) 


VW Cha 


8.1 (0.5) 


8.4 (0.3) 


115.0 (7.0) 


VZCha 


7.7 (0.6) 


7.7 (0.9) 


79.2 (2.4) 



° Value in parentheses indicates the spread in Piq for different J. 
Low-y includes selected lines up to P(14). High-y includes lines from 
P(22) to P(32). 



4.2. Excitation temperatures 
4.2.1. Rotational temperatures 

The extracted fluxes of the '^CO v = 1 - and '■'CO lines are 
used to produce excitation diagrams and thus calculate the gas 



rotational temperatures. This calculation is done by using the 
Boltzmann distribution, which assumes that the excitation in the 
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Fig. 14. Velocity shifts Vgas - Vstar as functions of Vstar- The gas velocity refers to the central velocity of the narrow component of 
the CO V = 1 - lines. 



gas can be characterized by a single excitation temperature and 
that the emission is optically thin; 



gi 



(1) 



and where the column density A^, of the upper level ; is de- 
scribed by, 



QAijhv 



(2) 



where the solid angle Q = D^^jcfi includes the emitting area 
Dem and the distance d to the object. The other quantities are 
the Einstein coefficient Ay [s '], the frequency v [cm '] and the 
measured line flux F-, [W m"^], the statistical weight gi of the up- 
per level, the total column density A^v for vibrational level v, the 
excitation temperature Tex and the rotational partition function 
QriTex)- Plotting the upper energy levels £, against fn{Fi/AijVgi) 
will give the rotational temperature Tiot by measuring the slope 
of the fitted line given by -1/rrot. 

The rotation diagrams for the '^CO lines are not linear, 
which suggests that the lines are optically thick. In addition, 
UV-pu mping can cause non-U nearity between the different tran- 
sitions (i Krotkovet allUQSOh . To overcome this problem, rota- 
tional temperatures have been estimated using only the '^'CO 
emission lines for the 4 sources, AS 205 A, DR Tau, RU Lup 
and S CrA A, for which sufficiently high S /N emission lines are 
detected, see Fig. [TS] and Table |7] The relative error of the line 
fluxes is taken to be 15% based on the noise in the spectra and 
uncertainties in the baseline fits. Only lines with J„ >1 are taken 
into account in the fit. 

While '^CO is significantly less abundant than '^CO, '^^CO 
lines may also be optically thick. '^CO lines with r ^1 have, for 
example, been seen in the disk around Herbig star AB Aur by 
[Blake & Boogert (2004). The inferred rotation temperatures of 
the 4 sources span 300-600 K, w hich is lower than 1 100-1300 K 
found for the T Tauri sample o f Naiita et al.l (l2003h but wi t hin th e 
temperature interval of 250-800 K given in lSalvk et all (l2009h . 
The discrepancy with Naiita et al. (2003) may arise from differ- 
ences in methodology. Naiita et al. (2003) fit '^CO data covering 
lines with upper energies levels of up to 10"* K, which they as- 
sume to be optically thin. They performed a linear fitting based 



on this assumption, but curvature in their rotation diagrams indi- 
cates that their '^CO lines are likely also optically thick, leading 
to uncertainties in their extracted temperatures. Another expla- 
nation for the temperature difference may be the lack of higher 
7- transitions in our data set. However, AS 205 A includes '^CO 
transitions up to £, - 4400 K, and still has a lower tempera- 
tureof 550 + 40 K. The temperature estimates by Salyk et al] 
(l2009l) based on '^CO should be more reliable since they took 
optical depth effects into account. In addition, a rotational tem- 
perature of 770 ± 140 K has been estimated for DR Tau using 
the C'^O lines, which is close to the temperature of 510 + 40 
K given by the '^CO lines (Fig. [TSll. If the two highest-7 '^CO 
lines are excluded from the fit, the '^CO rotational temperature 
increases to 610 ± 70 K, which is consistent with the C"*0 tem- 
perature within the error bars. These two high-7 '^CO lines are 
very weak and their fluxes may have been underestimated. 

We determine the '^CO/C'**0 ratio for DR Tau to examine 
whether the '^CO emission is indeed optically thin. The ratio 
is calculated by comparing the fluxes in the R(5) and P(16) 
Unes for '^CO with the R(5) and P(17) fines for C"*0, cho- 
sen because they are not blended and have similar J. This gives 
'^CO/C'*0 = 6.3 +1.0 which is close to the overall abundance 
ratio of '^CO/C'^O = 8.5 observed in the solar neighborhood 
d Wilson & Roodll 19941) . The optical depth of '^CO is estimated 
to be T = 0.3 + 0.2. Inferred column densities and emitting areas 
will be presented in Brown et al. (in prep.), based on combined 
fits to the I^CO and '^CO lines. 



4.2.2. Vibrational temperatures 

Vibrational temperatures are determined for the three sources, 
AS 205 A, DR Tau and S CrA A, for which high S /N '^CO 
data exist for several vibrational transitions (see Fig. [T6b . The 
vibrational temperatures are calculated using the same relation 
between the flux of the lines and the upper energy levels as 
presented in equations [1] and |2] The assumption made here is 
that every vibrational level has the same rotational temperature. 
The only difference is that in(Fi/AijV) instead of {n(Fi/AijVgi), 
is plotted vs the upper energy levels £, to extract the vibrational 
temperature Tyn, since the statistical weights of the vibrational 
levels are the same. One clear line, close in /-level and with- 
out strong line overlap, is chosen per ro-vibrational level (CO 
V = 1 - 0, V = 2 - 1 etc.) to extract the line fluxes. 
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Fig. 15. Rotational diagrams for T-Tauri stars AS 205 A, DR Tau, RU Lup and S CrA A. Blue stars represent '^CO lines and green 
squares C'^O lines. 

Table 7. Rotational temperatures Tiot derived from the '^CO lines, together with the vibrational temperatures T^ib- 



Source 


^rot 


[K] 




AS 205 A 


550 


±40 


1740 ± 150 


DRTau 


510 


±40 


1680 ± 140 


RULup 


360 


± 20 




SCrAA 


420 


±30 


1730 ± 140 




Fig. 16. Vibrational diagrams for AS 205 A, DR Tau and S CrA A. Triangles represents upper limits. 



The inferred vibrational temperatures are 1740 ± 150 K for 
AS 205 A, 1680 ± 140 K for DR Tau and 1730 + 140 K for 
S CrA A. The same fits have been done without including the 
'^CO V = 1 - lines since the lines within this band can be 
more optically thick than those for the higher vibrational tran- 
sitions. However, excluding the '^CO v = 1 - lines does not 
significantly change the results within the uncertainties. The up- 
per limits of the v = 5 - 4 flux have been estimated using F - 
1.065 Ipeak Ay, where the AV is taken to be the same as that of 
a CO V = 3 - 2 line and the amplitude Ipcait = 3cr/fo/n"'^^, where 
the parameter bin is the number of pixels within AV and cr is the 
rms of a clean part of the spectrum where the line of interest is 
expected. The detection level is set to 3cr. 

The vibrational temperatures of ~1700 K are in general 
higher than the rotational temperatures derived from the '^CO 
lin es of 300-600 K. Th e inferred temperatures can be compared 
to iBrittain et al.l (l2007h who find a lower rotational excitation 
temperature of ~200 K compared with a vibrational temperature 
of 5600 + 800 K for the disk ai'ound the Herbig Ae star HD 
141569. Similarly, van der Plas et al. (subm.) detect lower rota- 
tional temperatures of ~1000 K relative to vibrational values of 
-6000 - 9000 K for three other Herbig Ae/Be stars. Both pa- 



pers interpret the higher vibrational temperatures as caused by 
UV-pumping into the higher vibrational levels. 

The results presented here indicate that UV-pumping may 
also be an important process for populating the higher vibra- 
tional levels in disks around T Tauri stars, especially for the 
broad-based single peaked sources. The rotation/vibration tem- 
perature differences in Brittain et al. (2007) and van der Plas et 
al. (subm.) are, however, much larger than those derived here. 
This may be due to the stronger UV fields from Herbig Ae/Be 
stars than T Tauri stars in the wavelength range where CO is 
UV-pumped. 

Additional UV flux above the stellar photosphere can be pro- 
duced by accretion (IValenti et alil2000l) . Table[3]includes accre- 
tion luminosities taken from the literature for our T Tauri stars. 
These accretion luminosities are based on measurements of the 
hydrogen continua or Ho' line emission from the sources. The av- 
erage accretion luminosity for those broad-based single peaked 
lines for which accretion luminosities have been measured, is 
0.5 Lq. This value is higher than the average accretion luminos- 
ity of 0.1 Lq which is calculated for the rest of the sample of 
disks that are shown in Table [3] The strong Pfund-y6 lines ob- 
served in our data provide an additional confirmation of ongo- 
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ing strong accretion (see Fig. Higher accretion rates lead to 
higher UV fluxes, which, in turn, increase the UV pumping rate 
of CO into higher vibrational states. Quantitative estimates of the 
effect have been made by Brown et al. (in prep.) using the actual 
observed UV spectra for stars with spectral types A-K in an UV 
excitation model. Even if the enhanced UV from accretion for 
the T Tauri stars is included, however, the resulting vibrational 
excitation temperatures for a K-type star with additional UV due 
to accretion are still lower than those using an A-star spectrum. 
For the v = 1 - emission, the question remains how much 
UV-fluorescence contributes relative to thermal excitation. 

4.3. Lack of extended emission 



the central velocity of the line, however this can be explained by 
contamination from the companion star. The small feature seen 
for S CrA A may be a detection of extended line emission. VW 
Cha and VZ Cha have higher upper limits of the projected radial 
line extent of 10.4 and 6.3 AU respectively which is caused by a 
higher noise in the second moment of the continuum emission. 

Since the position angle of the slit is arbitrary with respect 
to that of the disks in our observations, an additional check on 
the lack of radial extended emission is provided by the different 
spatial line profiles of the stacked CO v = 1 - lines of AS 
205 A, DR Tau and S CrA A taken at 3 different rotational an- 
gles. All three show exactly the same profile. This means that the 
CO emission originates from within a few AU at many different 
angles around the sources. 
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Fig. 17. Top: The spectral line CO v = 1 - P(8) for AS 205. 
Bottom: The FWHM of the spatial profile of the above spectrum 
at each individual wavelength point. No increase in FWHM is 
seen at the position of the line. 

Information about the spatial extent of the CO emission from 
broad-based single peaked sources is an important ingredient for 
constraining models of their origin. The variance or formal sec- 
ond moment jt/2 of the spectral trace is calculated at each wave- 
length from the 2-dimensional spectrum, according to the fol- 
lowing equation. 

IJ2 = - Cf X (3) 

where jc, is the spatial position, is the flux at that point and C 
is the centroid position computed from C - 2x, x Fj/'EFj. 

For continuum emission, the average value of the second mo- 
ment reflects the spatial resolution provided by the AO system, 
assuming that the continuum is not spatially extended. If the line 
emission is spatially extented, the values of jj.2 will be larger 
than the noise in the continuum /i2 signal at wavelengths where 
the line emits. There is no significant detection of extended line 
emission for any of our sources, see Fig. [17] and Fig. [18] An 
upper limit on the radial extent of the line emission, Miune, is 
therefore calculated using the formula ARii„g^ < (A/?,„j„- + o-/^)^ 
-ARi„sir^, where AR^^^^^^ is equal to the mean of /i2/2 of the con- 
tinuum and (TAfi is the standard deviation of (yU2/2)'^^. The val- 
ues listed in Table [8]include an additional correction factor (1 -H 
C/L), where C/Lis the continuum/line flux ratio, because the to- 
tal centroid is diluted by the continuum flux. These results show 
that most of the projected CO emission must originate from 
within a few AU for AS 205 A, DR Tau, RU Lup and S CrA 
A. S CrA B and VV CrA A (not shown) have a feature close to 



5. Discussion 

Three possibilities for the origin of broad-based single peaked 
profiles are discussed below: a Keplerian rotating disk, a disk 
wind or a funnel flow. 

5.1. Rotating disk 

The results presented in ^3.21 show that broad-based sin- 
gle peaked lines (category C) cannot be reproduced using 
a Keplerian model with a power-law temperature profile. 
Increasing the local line broadening parameter AVmod to ~8 km 
s ' removes the double peak, but a larger local line broadening 
is not sufficient to give a good fit between the disk model and the 
data. As seen in § 13. 2. 31 a Keplerian disk model can fit the data by 
solving for a temperature structure. The best fits have a plausible 
temperature structure but require extended emission, possibly in 
combination with an increase in the local broadening parameter. 

The results of § l4.3l show that the bulk of the emission comes 
from within a few AU which means that a Keplerian model with 
a non-standard temperature profile including extended emission 
out to 30 AU cannot explain the broad-based single peaked line 
profiles. Even the model with enhanced turbulence and emission 
out to 10 AU becomes marginal, as illustrated by the expected 
second moments in these non-standard cases (Fig. [19]). The ver- 
tical error bars in Fig. [T9]represent the 3cr noise limit in AU for 
AS 205 A. Similar results hold for most of the other sources, 
although the Keplerian models with a non-standard temperature 
profile including an enhanced turbulence cannot be fully ruled 
out for sources like VW Cha and VZ Cha for which the limits 
on the spatial e xtent are larger. 

ICarr et al ] (12004 invoked a larger local line broadening of 
7-15 km s ' in their disk models to obtain a better fit to CO v = 
2-0 band head profiles observed toward various young stellar 
objects. As discussed in their paper, such values are much larger 
than the sound speed of a few km s ' in the inner disk, whereas 
models of turbulence in disks typically give velocities that are 
less than the sound speed (e.g.. lKlahr & Bodenheime3. l2006l) . In 
addition, these turbulent values are much highe r than those in- 
ferred from observations of gas in the outer disk (lOi etalll2006l: 
iHughes etalL 12009'). although the turbulence in a warm surface 
layer near the star may differ from that in the outer regions. A 
m ore extended dis cussion of turbulence in inner disks is given 
in ICarr et afl (12004 ). Even though we cannot fully exclude an 
increased level of turbulence for some of our sources, this is 
unlikely to be the sole explanation for the broad-based single 
peaked line profiles observed here. 

In summary, the above discussion shows that a pure 
Keplerian disk model does not fit the data. However, the es- 
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Table 8. The upper limit on the projected radial extent A/?/,,,^ of the line emission, together with AT?,, 



Source 




ARn„, [AU] 


AS 205 A 


15.1 : 


tO.l 


<2.1 


DRTau 


14.5 : 


tO.l 


<1.7 


RULup 


17.9 : 


tO.2 


<2.5 


S CrA A 


15.9: 


tO.2 


<2.0 


S CrA 








VV CrA~A* 


16.6: 


tO.5 




VW Cha 


23.4 : 


t2.3 


<10.4 


VZ Ciia 


20.7 : 


t 1.0 


<6.3 



" See text for definition 

* Contamination by companion star. 



timated rotational temperatures of 300 to 800 K for the broad 
centrally peaked sources are consistent with the temperature ex- 
pected for the inner parts of a disk. The symmetry of the lines 
and the small radial velocity shift between the gas and the star (a 
maximum of ~5 km s"') are additional indications that the emis- 
sion may originate in a disk. Therefore, it is likely that at least 
part of the emission originates from the disk but some additional 
process must contribute to these line profiles as well. 

5.2. Disk wind 

5.2.1 . Thermally launched winds 

Strong outflows are excluded due to the small velocity shifts of 
the molecular gas relative to the central star. Thermally launched 
winds are driven by irradiation of either FUV (far- ultraviolet), 
EUV (extreme ultraviolet ) , or X-ray photons ([Alexanden 
l200i iGorti & HollenbachL 120091: lErcolano & Owenl Isubm " 



Irradiation of the upper layers of the disk heats the gas to such 
high temperatures that the thermal energy exceeds the gravi- 
tational potential, leading to a flow off' the disk surface. The 
launching radius is roughly given by Rg = GM-^/V^^^, where Vorb 
is the orbital speed of the gas. The wind retains the Keplerian ro- 
tational velocity of the launch radius. 

If the disk is illuminated by EUV radiation, an HII region 
forms in the upper layer of the disk atmosphere. The ionized gas 
reaches a temperature of lO'* K, which corresponds to V ~10 km 
s"^ For a 1 Mq star, the launching radius is then Rg~ 5-10 
AU, and the wind is launched with a velocity of ~ 5 - 10 km s"' 
(lAlexandeil |2008). This small velocity shift is consistent with 
the small velocity shifts we see in the CO gas emission relative 
to the stellar radial velocity. However, an EUV-driven wind con- 
sists mostly of ionized gas, so that CO may be photodissociated. 
Detailed models of the CO chemistry are needed to test whether 
enough CO can survive in such a wind to produce the observed 
CO line profiles. 

In the X-ray and F UV photoevaporation models of 
iGorti & H ollenbach (2009), the evaporating gas is cooler than 
the EUV-iiTadiated gas. The expected launching radii of 10-50 
AU for X-ray photoevaporation and >50 AU for FUV photoe- 
vaporation are inconsistent with the lack of spatially-extended 
emission in our CO spectra. However, in models of photoevapo- 
ration by soft X-ray emission by lErcolano & O wen ( subr nj, the 
launching radius is closer to the star. lErcolano & Owen flsubm.h 
calculate profiles for atomic and ionized species within this 
model and predict symmetric, broad-based, single-peaked line 
profiles with small blueshifts of 0-5 km s"' . 

As for the EUV models, a detailed CO chemistry needs to 
be coupled with the wind model to test this scenario, but the fact 



that the wind is mostly neutral will help in maintaining some CO 
in the flow. 

In summary, EUV and soft-X-ray photoevaporation of the 
disk might produce some CO emission at low velocities. More 
detailed studies of the CO chemistry and velocity fields in the 
thermally-launched wind are needed to determine where CO can 
survive in the wind and the resulting emission line profile. 

5.2.2. Magnetically launched winds 

Magnetic fields of Classical T Tauri stars exert torques that 
are thought to be st rong enough to launch powerful winds. 
lEdwards et al.l (120061) describe two types of winds that are seen 
in absorption in the He 1 /110830 line, one with velocities of 
around 250 km s"' and another with lower speeds up to 50 km 
s"' . These velocity shifts are also commonly observed in optical 
forbidden lines, such as [O I] 6300 A. The high velocities are 
not consistent with the velocities seen in the broad-based single 
peaked CO emission, making this scenario highly unlikely. 

5.3. Funnel flow 

Accretion from the disk onto the star is thought to occur in a fun 



nel flow alon g strong dipol a r field lines (e.g., Hartmann et al 
1994; Bouvi eretal L l2007t lYang et al.L [2007; Gregory et si 



2008)). iNaiita et alj d200 3) rule out the possibility that the fun 
nel flow could produce some CO emission because the temper- 
atures in the fu nnel flow are expected to be ~ 3000 - 6000 K 
(iMartinl Il997 f). higher than the measured temperatures of the 
CO-emitting gas. Models of observed hydrogen emissi on fine 
profil es are consistent with th ese temper atures ([M uzeroll e et al.L 
200 it iKurosawa etall l2006l) . However, iBary et al.. {200§) sug- 
gest that the relative hydrogen emission fine fluxes may also be 
consistent with cooler temperatures of ~ 1000 - 1500 K, and 
cooler temperatures may be pre sent where the funnel flow con- 
nects with the disk. lNaiita et alJ (f2003) show that a modeled CO 
emission line profile from a funnel flow would have a double 
peaked line profile with no velocity shift (see Fig. 10b in their 
paper), which disagrees with the presence of the narrow emis- 
sion peak seen in the line profiles discussed here. Thus, funnel 
flows are also not likely to be the origin of the observed emis- 



sion. 



6. Conclusions 

Using CRIRES spectroscopy of CO (Av - 1) ro-vibrational 
emission lines around 4.7 //m from a sample of 50 T Tauri 
stars, we find that the line profiles can be divided into three 
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Fig. 18. The stacked square root of the formal second moment (jui)'^^ of the spectral trace for the '^CO emission lines. About 5-10 
clean Unes for each source were used for stacking. The square root of the second moment (pii)^^^ is multipUed with * = 0.2 for 
sources with a larger uncertainty in the second moment. The vertical shift for respective source is: AS 205 A = AU, DR Tau = 
-3.4 AU, RU Lup = -7.0 AU, S CrA A = -7.2 AU. For VW Cha and VZ Cha the shifts are 5.5 and 3.7 AU after multiplication with 
0.2. 



basic categories: A) narrow single peaked profiles, B), double 
peaked profiles and C) single peaked profiles with a narrow peak 
(~ 10 - 20kms"'), but broad wings, sometimes extending out 
to lOOkms"'. The broad-based single peaked line profiles are 
rather common, since they were detected in 8 sources in a sam- 
ple of about 50 T Tauri stars. These 8 sources have preferentially 
high accretion rates, show detections of higher vibrational emis- 
sion lines (up to v = 4) and have CO lines with unusually-high 
line-to-continuum ratios relative to all other sources within our 
sample. For at least one of the disks (S CrA B), the narrow com- 
ponent decreases faster for higher J Unes than the broad compo- 
nent, implying that the narrow component is colder. Generally, 
however, the two component fits give similar temperatures and 
central velocities. 

The broad-based single peaked emission originates within 
a few AU of the star with gas temperatures of ~300 - 800 K. 
The line profiles are symmetric and the gas radial velocity is 
close to the stellar radial velocity. These characteristics point to- 



ward an origin in a disk. However, unhke the CO profiles from 
other objects that are double-peaked or have only narrow peak, 
the broad-based single peaked line profiles could not be well fit 
using a Keplerian model with a power-law temperature profile. 
The fits are improved if a large turbulent width of ~ 8 km s"' 
is invoked, but the overall profile fit is still poor. Models with 
an unusual temperature distribution, perhaps with enhanced tur- 
bulence, provide a much better fit to the line profiles but require 
emission out to larger distances than observed. Thus, the hypoth- 
esis that this emission originates in a pure Keplerian disk needs 
to be questioned and an additional physical component needs to 
be considered. 

Several other scenarios are also ruled out. FUV radiation- 
driven winds have a launching radius of > 50 AU that is incon- 
sistent with the lack of spatially-extended emission in our data. 
Magneto-centrifrugal winds are observed to have blueshifts of 
50 - 250 km s"' , which are inconsistent with the lack of a signif- 
icant velocity shift in the symmetric CO emission lines. A funnel 
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Fig. 19. The second moment of the modeled spatial line profile for AS 205 (10 AU - O'.'OS at a distance of 125 pc). The model in 
the left plot includes extended emission out to 30 AU and Aymod=3 km s"' . The model in the right plot includes emission out to 10 
AU and an increased Ayn,od=7 km s"' (see Fig.|5]l. Error bars represent the 3cr noise limit for AS 205 A in Fig. [18] Both models use 
the observed position angle of 1 15° which is 50° off from the position angle of 165° of the AS 205 A disk tAndrews et al.t,2009.) . 
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flow is also unlikely because emission near the inner rim of the 
disk should have a double-peaked line profile, and because the 
temperatures in the funnel flow are expected to be higher than 
the calculated CO rotational temperatures. 

The most plausible explanation for the broad-based single- 
peaked line profiles is therefore some combination of emis- 
sion from the warm surface layers of the inner disk, contribut- 
ing to the broad component through Keplerian rotation, and a 
disk wind, responsible for the narrow component. A thermally 
launched disk wind, perhaps driven by EUV radiation or soft 
X-rays will have a small velocity shift of AV of 0-10 km s 
which is marginally consistent with the maximum detected ve- 
locity shift of AV ~5 km s"'. 

More information, including independent estimates of incli- 
nation angles, detailed models of CO chemistry and line profiles 
that originates in disk winds and the inner regions of disks, and 
spectro-astrometry of emission from these sources, is needed to 
distnguish the contributions from a disk (perhaps with enhanced 
turbulence) and a disk wind. These line profiles will be further 
discussed by Pontoppidan, Blake & Smette (subm.) using a com- 
bined disk and slow molecular disk wind model to analyze the 
spectro-astrometric data of these sources. 
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